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ABSTRACT 

Using the OSIRIS instrument installed on the 10.4-m Gran Telescopio Canarias (GTC) 
we acquired multi-color transit photometry of four small (R p 5 R$) short-period 
(P ^ 6 days) planet candidates recently identified by the Kepler space mission. These 
' observations are part of a program to constrain the false positive rate for small, short- 

PjJ . period Kepler planet candidates. Since planetary transits should be largely achro- 

matic when observed at different wavelengths (excluding the small color changes due 
to stellar limb darkening), we use the observed transit color to identify candidates 
as either false positives (e.g., a blend with a stellar eclipsing binary either in the 
background/foreground or bound to the target star) or validated planets. Our results 
include the identification of KOI 225.01 and KOI 1187.01 as false positives and the 
tentative validation of KOI 420.01 and KOI 526.01 as planets. The probability of iden- 
tifying two false positives out of a sample of four targets is less than 1%, assuming 
an overall false positive rate for Kepler planet candidates of 10% (as estimated by 
iMorton fc Johnso"n]|201l[h Therefore, these results suggest a higher false positive rate 
for the small, short-period Kepler planet candidates than has been theoretically pre- 
dicted by other studies which consider the Kepler planet candidate sample as a whole. 
Furthermore, our results are c onsistent with a recen t Doppler study of short-period gi- 
ant Kepler planet candidates (|Santerne et alj|2012j) . We also investigate how the false 
positive rate for our sample varies with different planetary and stellar properties. Our 
results suggest that the false positive rate varies significantly with orbital period and is 
largest at the shortest orbital periods (P < 3 days), where there is a corresponding rise 
in the number of detached eclipsing binary stars (i.e., systems that can easily mimic 
planetary transits) that have been discovered by Kepler. However, we do not find 
significant correlations between the false positive rate and other planetary or stellar 
properties. Our sample size is not yet large enough to determine if orbital period plays 
the largest role in determining the false positive rate, but we discuss plans for future 
observations of additional Kepler candidates and compare our program focusing on 
relatively faint Kepler targets from the GTC with follow-up of Kepler targets that 
has been done with w arm- S 'pit zer. 
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1 INTRODUCTION 

The Kepler space mission has discovered, to date, 61 tran- 
siting planets as well as ov er 2,000 planet candidates and 
2,000 eclipsing binary star s ( Batalha et al.ll2012i ; IPrsa et al.l 
l201ll : ISlawson et"afl 120111 ) Fl With such a vast number of 
planet candidates, it can be difficult to decide which to fo- 
cus follow-up efforts on. Recent studies have tried to address 
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this issue by estimating the false positive rate for the Kepler 
sample. Base d on the list of 1,235 K epler planet candidates 
published by iBorucki et all (|2011bh . it has been predicted 
that as many as ^95% of these candidates are true planets 
(|Morton fc Johnson! [20TI . However, previous studies have 
not taken into account how different subsets of Kepler tar- 
gets may have different false positive rates. For example, 
there is a rapid rise in the number of detached eclipsing bi- 
nary stars that have been discovered by Kepler at orbital pe- 
riods of less tha n ~ 3 days, and s u ch systems can mimi c plan- 
etary transits (|Prsa et al.ll201ll ; ISlawson et al.l 120111 ). This 
suggests that there may be corresponding changes in the 
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false positive rate with orbital period. Because the probabil- 
ity of observing a transit event increases as the orbital period 
of the planet decreases, many of Kepler's planet candidates 
have short periods. Thus it is necessary to be cautious when 
estimating false positive rates over the whole Kepler sam- 
ple. While observational studies with waxm-Spitzer support 
predicti ons of low false pos itive rates over the entire Kepler 
sample (iDesert et alj|2012h . biases in target selection can af- 
fect observationally-constrained false positive rates (see 33 
for further discussion). We also note that a recent imag- 
ing study has found that nearly 42% of their sample of 
98 Kepler planet candidate hosts has a visual or bound 
companion within 6 arcseconds of the target star (Lillo- 
Box et al., in preparation). Such studies emphasize the need 
for follow-up imaging to exclude blend scenarios imitating 
planet candidates or contaminated transit depths. 

The false positive scenarios we consider in this paper 
are those that result from stellar eclipsing binaries that are 
either in the background (or, in rare cases, foreground) or 
bound to the target star and are not always easy to identify 
with Kepler due to the flux from the different stars being 
blended together within Keple r's photometri c apert ure (~6 
arcsec). As discussed by, e.g., I Colon fc Ford! (^Olll ). differ- 
ent techniques can be used to eliminate many, but not all, 
blends. Multi-color transit photometry is an efficient method 
for recognizing blends that cannot be spatially resolved, as 
measuring the transit depth in different bandpasses (i.e. the 
transit color) allows one to test the planet hypothesis. This 
is possible since the magnitude of the transit color changes 
as long as the blended stars have significantly different col- 
ors. 

IColon fc Fordl (|201 ll) presented multi-color transit pho- 
tometry of a Kepler target , KOI 565.01, that was first 
announced by iBorucki et alj (l2011aT l to be a super-Earth- 
size planet candidate but later recognized as a likely false 
positive due to measurements of a centroid shift away 
fro m the location of the tar get on the CCD du ring tran- 
sit l|Borucki et al.ll2011bl) . In lColon fc Fordl (|201ll ). we used 
near-simultaneous multi-color observations acquired using 
the narrow-band tunable filter imaging mode on the Opti- 
cal System for Imaging and low Resolution Integrated Spec- 
troscopy (OSIRIS) installed on the 10.4-m Gran Telescopio 
Canarias (GTC) to confirm that KOI 565.01 is indeed a 
false positive, as we both resolved a stellar eclipsing binary 
~15 arcsec from the target and measured a color change 
in the "unreso l ved" target+eclipsing binary system. Thus, 
IColon fc Fordl d 20111 ) demonstrated the capability of the 
GTC/OSIRIS for efficient vetting of planet candidates via 
its capabilities for near-simultaneous multi-color photome- 
try within a single transit event. 

In this paper, we present observations of four Kepler 
planet candidates specifically selected to have small radii 
and short orbital periods (R p ^ 5 ii® and P ^ 6 days). 
Measuring the false positive rate for this extreme subset of 
planet candidates will allow us to test if there is a correlation 
between the false positi ve rate and dif f erent planetary and 
stellar properties. As in Co lon fc Fordl l|201ll ). the observa- 
tions presented here were acquired with the GTC/OSIRIS. 
However, these observations used broadband filters in lieu 
of the narrow-band tunable filters in order to collect more 
photons and to obtain greater wavelength coverage (and 
thereby probe greater color differences) . In ij2] we discuss 



our target selection criteria, and we discuss the correspond- 
ing observations for our four targets in Sj3] In Sj4]and fj5]wc 
describe our data reduction and light curve analysis pro- 
cedures, and we present results for each target in i|6] We 
include a discussion of our results and how they relate to 
the distribution of eclipsing binaries that have been discov- 
ered by Kepler as well as previous estimates of the false 
positive rate for the Kepler sam ple in jjTl In particu l ar, we 
discuss theoretical estimates from lMorton fc Johnson! (1201 if) 
and o bser vational constra i nts fr om studies bv lDesert et al.l 
(2012) and San terne et al.l (|2012T l. Finally, we summarize our 
results and conclusions in ijEl and we also discuss our plans 
for future observations of additional Kepler targets with the 
GTC. 



2 TARGET SELECTION 

Recent studies (e.g. jLissauer et~al ]|201ll . l2012l ) have demon- 
strated that a majority of the planet candidates in Kepler's 
multi-planet candidate systems should in fact be real plan- 
ets. Considering these studies, for our program we target 
on ly those candidates fou nd in single systems as determined 
bv lBorucki et al] (|2011bh . From the sample of candidates in 
single systems, we selected targets based on the following 
criteria: 

• orbital period (P) ^ 6 days 

• planet radius (R p ) ^ 6 -R® 

• transit depth (at center of transit; 8) > 500 ppm 

• transit duration (first to fourth contact; r) ^ 2.5 h 

• Kepler magnitude (Kp) < 15.5 

• vetting flag > 1 

As discussed in the previous section, there is a significant rise 
in the number of detached eclipsing binary stars compared 
to planet candidates at short orbital periods. This is also 
illustrated in Fig. [T] where we show histograms of the num- 
ber of planet candidates, detached eclipsing binaries, and 
all "other" eclipsing binaries (i.e. all bin aries that are not 
listed as detached in ISlawson et al.1 l201ll ) as a function of 
orbital period. Therefore, we focus on constraining the false 
positive rate only for planet candidates with short orbital 
periods, where the presence of eclipsing binaries is greatest. 
Our decision to focus on small candidates is due to there be- 
ing a dominant population of Neptune-size or smaller can- 
didates in the Kepler sample. We set constraints on the 
transit depth, transit duration and Kepler magnitude due 
to limitations of our observing technique, so as to measure 
the transit depth and color precisely. Furthermore, we want 
to do this while maintaining a reasonable observing cadence, 
hence the limits on the Kepler magnitude of our targets. Fi- 
nally, we exclude candidates that have vetting flags of 1, as 
those have been previously confirmed as planets. 

Along with the above constraints on the planetary and 
stellar properties, we have several observational constraints 
due to limitations of the GTC/OSIRIS. First, we rule out 
candidates that have no observable transits from the loca- 
tion of the GTC (which is located on La Palma at the Ob- 
servatorio del Roque de los Muchachos). Then, we exclude 
those that did not have multiple transits observable during 
bright or grey time (during the 2011 observing season). We 
also exclude transit events that occur outside altitudes of 
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Figure 1. Histograms illustrating the number of Keple r planet 
candidates (solid blue line; based on lBorucki et al .1201 lbh and de- 
tached and all "other" eclipsing b inaries (dotted red a nd dashed 
black lines, respectively; based on lSlawson et al1l201ll) as a func- 
tion of orbital period. We also show a histogram (dash-dot green 
line) based on the recently released list of 2,321 KOIs from 
iBatalha et al The righthand y-axis shows the percentage 

of systems relative to the total number of target stars observed 
by Kepler in Ql. Note that only systems with orbital periods of 
up to 10 days are shown. While there is a substantial number of 
"other" (primarily semi-detached and overcontact) eclipsing bina- 
ries at the shortest orbital periods (which may be less likely to be 
flagged as potential planets), there is still a greater-to-comparable 
number of short-period detached eclipsing binaries compared to 
short-period planet candidates, which indicates the potential for 
eclipsing binaries to infiltrate the short-period planet candidate 
sample. The orbital periods of the KOIs we observed are indicated 
with vertical long-dashed gray lines and are labeled accordingly. 



Figure 2. Radius versus orbital period of Kepler planet candi- 
dates. The locations of the four KOIs we observed are indicated 
with filled red squares and are labeled accord ingly. The locations 
of the targets observed with warm- Spitzer bv lDesert et al.l ||2012| ) 
are indicated with filled blue squares. The filled gray circles and 
th e open gray circles re spect ively represent all K OIs presented 
bv lBorucki et al] l|2011bri and lBatalha et al . (2012) out to orbital 
periods of 200 days. Note the prominence of KOIs smaller than 
~ 6 R®. 



ties of our targets in comparison to the sample of 1,235 KOIs 
and the corresponding 997 host stars from which our tar- 
gets wer e chosen. The samp le of KOIs observed with warm- 
Spitzer l|Desert et al.l [2013 ) is also illustrated in these fig- 
ures for comparison (we discuss the Spitzer sample in more 
detail in Sj7T}. Our targets and the corresponding observa- 
tions are described in detail in 33] below. 



~35-72 degrees, so as to avoid high airmass and vignetting 
that occurs at high altitudes due to limitations of the GTC 
dome. 

Finally, we did not check the Kepler light curved for 
each target that fit all of the above criteria prior to acquir- 
ing observations. However, as we discuss in f|8] for future 
observations we plan to check the Kepler light curves for 
secondary eclipses and V-shaped eclipses, either of which 
could indicate a priori that the target is not a transiting 
planet but instead a stellar eclipsing binary. 

These selection criteria led us to acquire observations 
of four planet candidates between April and September 
2011: KOI (Kepler Object of Interest) 225.01, 420.01, 526.01 
and 1187.01. In Table [T] we list so me of the properties o f 
each of our ta r gets as determined bv lBorucki et al] (l2011bl ). 
iBorucki et all (|2011bl) flagged KOI 225.01 as possibly hav- 
ing ellipsoidal variations. KOI 1187.01 has t he shortest pe- 
riod of all the planet candidates announced bv lBorucki et al.l 
l|2011bD . In Figures [3 E] and U we illustrate different proper- 



3 OBSERVATIONS 

We acquired photometry of each target and several nearby 
reference stars around the predicted time of a transit event. 
For each observation, we used the GTC/OSIRIS to ac- 
quire near-simultaneous multi-color photometry by alternat- 
ing between two broadband "order sorter" filterqj that were 
custom made for OSIRIS; 666T36 nm and 858T58 nm. Both 
bandpasses were specifically chosen so as to minimize effects 
of telluric absorption and emission while also allowing for 
ample wavelength coverage. In order to decrease dead time, 
all observations used lxl binning, a fast pixel readout rate 
of 500 kHz and a single window located on one CCD chip. 
The size of the window varied for each observation, but each 
window was large enough to contain the target and several 
reference stars. We note that all observations except those 
for KOI 1187 were conducted in queue (service) mode. In 
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Table 1. KOI Properties 



Star Planet Candidate 



KOI 


KIC 


Kp 


T e// (K) 


b (deg) 


P (days) 


r (hr) 


5 (ppm) 


Rp (R<s) 


Rpl R* 


225 


5801571 


14.784 


6037 


9.177 


0.838598 


1.2452 


2571 


4.9 


0.04932 


420 


8352537 


14.247 


4687 


16.412 


6.010401 


2.2582 


2700 


4.3 


0.0474 


526 


9157634 


14.427 


5467 


12.426 


2.104719 


1.7639 


926 


2.6 


0.0305 


1187 


3848972 


14.489 


5286 


10.846 


0.3705285 


0.7778 


1835 


2.5 


0.03961 



All values are from lBorucki et~al ] l|2011bh . The KIC number refers to the Kepler Input Catalog number for 
each target. Note that b is the Galactic latitude of the KOI host star, r is the transit duration and <5 is the 
transit depth as measured in the Kepler bandpass. 
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Figure 3. Signal-to-noise ratio (SNR) per transit as a function of 
orbital period for Kepler planet candidates. The SNR per transit 
for each KOI was computed from the transit depth, the transit 
duration and either the 6-hour combined differential photometric 
precision (CDPP) from Q3 data (for the Borucki et al. 2011b list) 
or the 6-hour CDPP from Q1-Q6 data (for the Batalha et al. 2012 
list). The colors and symbols are the same as in Figure [2] 

the following sections, we describe specific details regarding 
each target and its respective observations. 

3.1 KOI 225.01 

We observed the transit event of KOI 225.01 on 2011 April 
13 under clear conditions and during dark time, with ob- 
servations beginning at 03:10 UT and ending at 06:12 UT, 
during which time the airmass ranged from ~ 1.66 to 1.07. 
The exposure time was 50 s (for both bandpasses), with 
^40 s of dead time following each exposure. Seeing varied 
between 1.2 and 2.0 arcsec until about 05:00 UT, at which 
point the seeing improved to < 1 arcsec and a slight defo- 
cus was introduced in order to avoid saturation. A diffuse, 
dark band was present in all the images (a result of a very 
bright star located just outside the CCD), so the target was 
positioned on the CCD such that it was outside this region 



Figure 4. Kepler magnitude (Kp) versu s effective tem p erature 
for the 997 KOI host stars published in iBorucki et al.l J2011bll 
and th e additional 926 KOI host stars published in lBatalhaet al.l 
The colors and symbols are the same as in Figs. |21 and l3l 

and also avoided bad pixels in the center of the CCD chip. 
During the observations, the target's centroid coordinates 
shifted by < 9 pixels in the x-direction and ~ 2 pixels in the 
^-direction. Several images were lost due to technical issues, 
and a few of the images towards the end of the observations 
were discarded due to the beginning of twilight. 

3.2 KOI 420.01 

Observations of the 2011 September 13 transit of KOI 420.01 
took place from 21:48 UT (2011 September 13) to 01:30 UT 
(2011 September 14). Observations occurred during bright 
time and under photometric conditions, and the airmass 
ranged from 1.07 to 1.84. The exposure time was set to 40 s, 
with a corresponding 20 s of dead time following each expo- 
sure. The seeing roughly ranged from 1.1-1.8 arcsec through- 
out the observations, though at the beginning of the ob- 
servations the seeing conditions changed drastically enough 
from one image to another that some of the reference stars 
were saturated. We explicitly exclude any saturated refer- 
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ence stars from our analysis (see ij5]for further details). The 
seeing stabilized and improved after about 23:00 UT, enough 
so that a slight defocus was implemented to avoid satura- 
tion. The target's centroid coordinates shifted by £ 4 pixels 
in either direction during the observations. 



3.3 KOI 526.01 

The 2011 September 11 transit of KOI 526.01 was observed 
under photometric conditions and during bright time. Ob- 
servations began at 21:37 UT on 2011 September 11 and 
ended at 00:00 UT on 2011 September 12. The airmass 
ranged from 1.05 to 1.22 and the seeing was stable between 
0.6 and 0.8 arcsec. A slight defocus was implemented, yield- 
ing a defocused FWHM of 0.8-1.0 arcsec. An exposure time 
of 10 s was used, with 20 s of dead time between exposures. 
As in observations of KOI 225.01 f i|3.1[) . a dark band caused 
by a bright star outside the CCD window was present in all 
images, so we placed the target star appropriately far away 
from the band so that the photometry would not be af- 
fected. The centroid coordinates of the target shifted by less 
than 4 and 2 pixels in the x- and y-directions, respectively. 
Twice during the observations the primary mirror segments 
lost alignment and produced distorted images: first around 
22:30 UT and again at 00:00 UT. We discarded a few im- 
ages from the first instance, but the issue with the mirror 
could not be readily fixed during the second instance, and 
the observations were forced to end early, around the time 
of the transit egress. 



3.4 KOI 1187.01 

We observed the 2011 June 12 transit event of KOI 1187.01, 
with observations taking place from 22:31 UT (2011 June 11) 
to 02:01 UT (2011 June 12). Observations took place dur- 
ing bright time and under clear conditions, and the airmass 
ranged from 1.90 to 1.05. Seeing was excellent and ranged 
from 0.6 to 0.9 arcsec during the night. Due to the excellent 
seeing conditions, the exposure time was set to 5 s. There 
were 24 s of dead time following each exposure. During the 
observations, the target's centroid coordinates shifted by ^ 
4 pixels in either direction. 



4 DATA REDUCTION 

We used standard iraf procedures for bias subtraction and 
flat-field correction (using dome flats taken for each band- 
pass) for each target. Due to nonuniform illumination by 
the lamp, we added an illumination correction to the final 
flat-field image for each target (except KOI 225). The il- 
lumination correction was performed using the iraf task 
mkillumflat within the noao.imred.ccdred package, which 
removes the large scale illumination pattern from the flat 
field by smoothing the fiat field image. For KOI 225, the 
final flat-field for the 666 nm filter showed a strong gradient 
after performing an illumination correction, so we chose to 
use the fiat-field image as it was. To be consistent, we did 
not perform an illumination correction on the 858 nm final 
flat-field for KOI 225 either. We used the IDL Astronomy 



User's LibrarjQ implementation of daophot (|Stetsonlll98"7h 
to perform aperture photometry on each target and several 
nearby reference stars as well as on other potential sources 
of the transit signal (i.e. stars within ~ 20 arcsec of each tar- 
get). For each target, we tested different apertures and chose 
a final aperture based on that which resulted in the smallest 
scatter in the baseline (out-of-transit) flux ratios (i.e., the 
target star flux divided by the ensemble reference star flux). 
For KOI 225, 420, 526 and 1187 our final aperture was 25, 
23, 15 and 14 pixels (equivalent to approximately 3.2, 2.9, 
1.9 and 1.8 arcsec). Sky background subtraction took place 
during the aperture photometry process, with annuli chosen 
to be far enough away from each star that the flux from a 
given star would not be included within the sky annulus. 
While these procedures were performed for each bandpass 
separately, we considered the results for each bandpass and 
used the same aperture and sky annulus for each data set for 
a given target. Note that in many cases a smaller aperture 
and sky annulus were used for the stars within ~ 20 arc- 
sec of a given target due to the small separations between 
some of the stars and the potential for blending. For refer- 
ence, we present a reduced science image showing the field 
of view around KOI 1187 in Figure [5] Once aperture pho- 
tometry and sky subtraction were completed, we proceeded 
with light curve analysis as discussed in ij5] below. 



5 LIGHT CURVE ANALYSIS 

For each target, we computed light curves for each bandpass 
by dividing the total flux measured from the target by the to- 
tal weighted flux of an ensemble of reference stars, where the 
flux from each reference star was weighted accordingly after 
discarding outlying flux values (resulting from either tech- 
nical issues or variable sky conditions). During the analysis, 
several reference stars were found to be obviously variable 
or saturated, so we excluded those stars from the analysis. 
This resulted in us using 10, 4, 2 and 9 reference stars in the 
field near KOI 225, 420, 526 and 1187 to compute the ref- 
erence ensemble flux. Next, each light curve for each target 
was normalized to the mean baseline flux ratio as measured 
in each bandpass. The light curves were then corrected for 
changes in airmass as well as for drifts in the centroid co- 
ordinates of the target on the CCD and the sharpness of 
the target's profile [(2.35/FWHM) 2 ]. The latter corrections 
were done via ex t ernal parameter decorrelation (see, e.g., 
iBakos et ail 120071 . |2010|) . We note that due to the lack of 
data for KOI 526 during the transit egress and post-transit, 
our attempts to correct those light curves produced skewed 
results. Therefore, all further analysis for KOI 526 is based 
on the uncorrected light curves. 

The corresponding observation times for each light 
curve for each target were computed from the UTC times- 
tamps given in the image headers. We converted the 
UTC time at mid-exposure to Barycentric Julian Dates 
in Barycentric Dynamic al Time (BJD_ T DB) via an online 
calculator described in lEastman et~aH (|2010T ) in order to 
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Figure 5. Image of the field of view around KOI 1187 as ac- 
quired using the GTC/OSIRIS and the 666/36 nm order sorter 
filter. The target is labeled accordingly and is contained within 
a black circle equivalent to twice the aperture used in the reduc- 
tion process. The unlabeled stars contained in black circles are 
the nine reference stars used in our analysis. The larger red circle 
around KOI 1187 has a radius of ~ 20 arcsec, and the six brightest 
stars located within this circle (excluding the target) are the stars 
that we investigated as potential sources of the transit signal. 

matc h the time coordinat e system that the ephemerides 
from IBorucki et "all (|2011bl ) are given in0 

The photometric uncertainties were computed from the 
photon noise of the target star and the reference star en- 

5 The calculator is available at 

: / /astroutils. astronomy.ohio-state.edu/time/utc2bjd. html 



Table 2. Photometric Precisions 



KOI 


C666nm (ppm) 


0-858nm (PP m ) 


225 


476 


540 


420 


501 


428 


526 


1223 


1100 


1187 


1447 


1295 



semble, the noise in the sky background around the target 
and each reference used in the ensemble, and scintillation 
noise. The resulting median photometric uncertainties for 
each light curve for each target are given in Tabled In each 
case, the photon noise of the target is the primary source 
of error. Finally, in order to investigate the presence of red 
noise in our data, we computed the standard deviation of the 
flux ratios for different bin sizes and found that all our data 
follow the trend expected for white (Gaussian) noise. Despite 
this result, we also investigate potential residual systematics 
via a "prayer-bead" analysis, which we discuss below. 

After correcting the light curves for each target (except 
for KOI 526 as discussed above), we proceeded to fit syn- 
thetic mo dels to each light curve, following the approach 
taken by IColon et al.l (|2010l ). We assumed each observed 
transit event was a planetary transit , and we used th e plan - 
etary transit light curve models from lMandel fc Agoll (|2002l ) 
to fit limb-darkened models to our data. Specifically, for each 
target, we fit for the following parameters: 

• time of mid-transit (to) 

• transit duration (first to fourth contact; r) 

• impact parameter (a cos i/ J?*) 

• planet-star radius ratio (R p /R*) 

• two limb darkening coefficients (ci and C2) 

• baseline flux ratio 

• (linear) baseline slope 

The limb darkening coefficients that we fit for are defined as 
ci = Mi + U2 and C2 = Mi — U2, where Mi and M2 are linear 
and quadratic limb darkening coefficients. For KOI 526, we 
note that we held t he tran sit duration fixed to the value 
from IBorucki et all l|2011bl ) in order to fit a model to our 
partial transit. 

Initial guesses for the mid-transit time, transit dura- 
tion, impact parameter and radius ratio were based on the 
values and their cor responding uncertainties as given in 
IBorucki et ail (l2011bl ). Values for the limb darkening coeffi- 
cients were interpolated from the IClaret fc BloemenT (|201ll ) 
models for the Sloan r and z filters and are bas ed on the 
stellar parameters given in IBorucki et all (|2011bl ). Due to 
the potentially larger than estimated uncertainti es in the 
stellar parameters (see, e.g., IBorucki et al.ll2011bl ). we did 
not allow the limb darkening coefficients to be completely 
free parameters in the fitting process. Rather, we kept the 
coefficients fixed at self-consistent values during the fitting 
process, but we tested a range of fixed values for the coeffi- 
cients. Similarly, we tested a range of initial guesses for the 
other parameters based on the uncertainties for the plane- 
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tary and stellar parameters. Best-fitting models were iden- 
tified via a Levenberg-Marquardt minimization scheme0 

Our light curve fitting procedure for each target is as 
follows. First, we fit models to each light curve individually, 
corrected the data against the best-fit baseline flux ratio and 
slope, subtracted the best-fit models and discarded any data 
points lying greater than 3<r from the residuals. Then, we 
used the corrected light curves and fit models to them in a 
joint analysis, where we forced the impact parameter, transit 
duration, mid-transit time, baseline flux ratio and baseline 
slope to be the same for both light curves. However, we al- 
lowed different values for the radius ratio and limb darkening 
coefficients to be fitted to the different light curves. The re- 
sults from the joint analysis were then used to correct the 
individual light curves as in the first step described above. A 
final joint analysis was then applied to the "final" corrected 
light curves. During this final stage, we also perfo r med a 
"prayer-bead" analysis as was done in Col on et alJ £2010). 
Specifically, we performed a circular shift on the residuals 
for each light curve (computed after removing the best-fit 
model) and constructed synthetic light curves by adding 
the shifted residuals back to the best-fit model. The joint 
analysis described above was applied to each synthetic light 
curve, and we use the dispersion of the best-fit parameters 
to calculate uncertainties on each parameter. This accounts 
for any additional systematic noise sources in the data. We 
present results from our light curve analysis in Jj6j 



6 RESULTS 

We present the light curves and the corresponding best-fit 
models for each target in Figures [6l [8] and [9] While not 
shown here, we had also generated light curves for the poten- 
tial sources of each transit signal (i.e. stars that could have 
been blended with the target within Kepler's aperture) fol- 
lowing a similar procedure as described above. We generated 
light curves for fO, 3, 3 and 6 stars within ~ 20 arcsec of 
KOI 225, 420, 526 and 1 187. Upon visual inspection, we 
found that none of these light curves showed a transit sig- 
nal during the time of the transit event, indicating that the 
transits that we observed either occur due to an object tran- 
siting the target star or an object transiting an unresolved 
star that is blended with the target star. We also present 
the transit color (666 nm — 858 nm) of each target in the 
bottom panel of Figures [6r-[9l which was computed by tak- 
ing the average of each pair of flux ratios in the 666 nm light 
curve and dividing by the corresponding points in the 858 
nm light curve. Therefore, in Figures [IJHHl a positive color 
indicates a redder transit and a negative color indicates a 
bluer transit. 

In Table(3] we report best-fit model parameters for each 
target based on the model with the smallest \ 2 value, but we 
note that the parameter uncertainties are based on the full 
set of best-fit values and their corresponding uncertainties as 
found during the "prayer-bead" analysis. Specifically, for a 
given parameter, the upper error bar is a result of subtract- 
ing the best-fit value (from the model with the smallest \ 2 
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Figure 6. Light curves, residuals, and color based on observations 
of the 2011 April 13 transit of KOI 225.01 acquired at 666 nm (in 
blue) and 858 nm (in red). The circles are observations and the 
solid curves are the corresponding best-fit models (see text for 
further details). Representative photometric error bars for each 
light curve are shown on the left-hand side of the top panel. The 
black horizontal dashed line illustrates the depth of the transit as 
measured in the Kepler bandpass (see Table [TJ. Residuals from 
the fits for each light curve are also shown in the top panel, with a 
horizontal dotted line indicating the level at which the residuals 
were offset (for clarity). The bottom panel shows the color as 
computed from the two light curves as well as a representative 
error bar on the left-hand side of the panel and a horizontal dotted 
line that illustrates a color of zero (for reference). 
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6 We specifically used mp fit fun, which is publicly available at 
http://www.physics.wisc.edu/~craigm/idl/idl.html 



Figure 7. Same as Figurc[6] but for the transit of KOI 420.01 as 
observed on 2011 September 13. 
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Figure 8. Same as Figurel6l but for the transit of KOI 526.01 as 
observed on 2011 September 11. 
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Figure 9. Same as Figure[6] but for the transit of KOI 1187.01 
as observed on 2011 June 12. 



value) from the maximum sum of a fitted value and its asso- 
ciated measurement uncertainty as determined from the full 
set of models computed during the "prayer-bead" analysis, 
and likewise for the lower error bar. In general, we find that 
the formal la errors for the best-fit parameters are compara- 
ble to those errors computed based on the results from the 
"prayer-bead" analysis, which indicates that any residual 
systematics in the data (i.e. those not removed by external 
parameter decorrelation or airmass corrections) have a neg- 
ligible effect on our results. Furthermore, the distribution of 



each best-fit light curve parameter over all permutations is 
smaller than the uncertainties given in Table [3j which fur- 
ther indicates that the results presented here are robust. 

Finally, as suggested by the referee, we consider whether 
we can use astrometry from our images to provide additional 
constraints on the properties of the targets we observed. A 
preliminary astrometric analysis does not suggest that any 
of the KOIs we observed are due to a blend with another 
star that was resolved by our observations. Furthermore, 
despite the large aperture used by Kepler, the astrometric 
precisions from Kepler are extrem ely high (< . 004 a rcsec- 
onds for a single 30 min exposure; iMonet et al.ll201ol ). and 
to the best of our knowledge, no significant centroid shifts 
were measured by Kepler for any of our KOIs. 

We discuss specific results for each target individually 
in the following sections. 



6.1 KOI 225.01 

The light curves presented in Figure [6] are prominently V- 
shaped, which alone suggests a possible non-planetary tran- 
sit event. Considering the shape of the light curves and the 
significant difference in the transit depths as measured in 
the GTC bandpasses, we deduce that KOI 225.01 is most 
likely a stellar eclipsing binary system. Specifically, we find 
that the best-fit planet-star radius ratios differ at a level 
of > 11a, which is a clear indication that this object is 
not a planet. The transit depth in the Kepler bandpass is 
also significantly different from those measured in the GTC 
bandpasses. In particular, the transit depth is notably differ- 
ent between the Kepler bandpass and the blue GTC band- 
pass (around 666 nm) , even though they probe fairly similar 
wavelength regimes. We primarily attribute this difference 
to dilution in the Kepler photometry, since there is a rela- 
tively bright star near KOI 225 that appears to have slightly 
contaminated one of the pixels u s ed in t he Kepler photom- 
etryQ As noted in lBorucki et al.l (|2011bh . this target's light 
curve had possible ellipsoidal variations, also suggesting that 
the system contains an eclipsing binar y. Furthermor e , after 
our observations had been conducted, ISlawson et all (|201lD 
listed this targe t in th eir eclipsing binary cataloglf] Finally, 
lOfir fc Dreizlerl l|2012h recently conducted an independent 
analysis of the Kepler data set and identified KOI 225.01 as 
a false positive due to the detection of significant differences 
between the odd and even eclipse events. 

Since we measure a deeper eclipse in the bluer GTC 
bandpass, this indicates that during the eclipse more blue 
light is blocked, and therefore the secondary (eclipsing) com- 
ponent is redder than the primary star. In Table [3] we pro- 
vide an updated ephemeris and eclipse duration, though we 
note that likely due to the non-planetary transit light curve 
shape, the best-fit model yielded an impact parameter that 
reached the upper boundary limits during the fitting pro- 
cess. Therefore, the resulting best-fit light curve parameters 



7 As determined from target pixel files downloaded from 
http: / /archive, stsci.edu/kcplcr/ 

8 We refer the reader to Slawson et al. (2011) for further details 
on how some KOIs were rejected as planet candidates and subse- 
quently added to the eclipsing binary catalog. 
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Table 3. Best-Fit Model Parameters 



KOI 


to 


T 


a cos i/Ri, 


Rp/ Ri, 


Rp/ R+ 


Cl 


C2 


Cl 


C2 




(BJD-2454900) 
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(858 nm) 
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n iQia+0.0026 
U.1318_o.ooi5 
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nfi^^+O 0035 
U.UOdO_ 00 25 

0924+ - 0174 
u.uyz^_ 0173 


0.6548 
0.7581 
0.6920 
0.6895 


0.0152 
0.4729 
0.2466 
0.2387 


0.5034 
0.5837 
0.5339 
0.5370 


-0.0398 
-0.2011 
0.0869 
0.0886 



Note that the time of mid-transit, to, is technically given in BJD_TDB (Barycentric Julian Date in Barycentric Dynamical Time). See text 
for further details. 



are likely somewhat skewed, but this does not change the 
signifi cance of our r esults . 

ISlawson et alj (|201ll) report a pe ri od tha t is twice as 
long as that fo und by IBorucki et al.1 |2011bl ). If correct, 
then either what IBorucki et al. ( 2011bl ) believed were tran- 
sit events were actually a combination of primary and sec- 
ondary eclipses or the eclipsing binary has a large enough 
inclination that all eclipse events are primary eclipses and 
the secondaries are unobservable fr om our line-of - sight . If 
we assume that the period found bv lSlawson et al.1 (|201ll ) is 
correct (and was determined, for instance, by measuring dif- 
ferent depths for successive transits, which led to this KOI 
being rejected), then it is necessary to consider how this af- 
fects any correlation between the orbital period and the false 
positive rate for Kepler targets. We refer the reader to SJ7] 
for further discussion. 



6.2 KOI 420.01 

As illustrated in Figure [7] the transit depths as observed 
in the GTC bandpasses are comparable, and there is not 
a significant change in the measured color during transit. 
The planet-star radius ratios measured from the GTC light 
curves are consistent within ~ 2.8a, and we did not resolve 
any potential sources of the transit signal within ~ 20 arc- 
sec of the target, so we deduce that this target is not a false 
positive and is instead a validated planet 3 The best-fit pa- 
rameters from our light curve models are given in Table 
and include an updated ephemeris, transit duration, and im- 
pa ct parameter. We not e that the impact parameter given 
in IBorucki et al.l ll2011bl) h as an associated uncertainty of 
1, and Batalha et al. ( 20121 ) report an impact parameter of 
0.57 with a relatively large uncertainty of 0.52, so we provide 
a much stronger constraint on this parameter here. Finally, 
based on our measured planet-star radius ratios and the re- 
vised stellar ra dius (0.69 Rq compar ed to th e value of 0.83 
Rp, re ported bv lBorucki et al.ll2011bl ) given in lBatalha et al.l 
|2012T ). we find KOI 420.01 to have a radius that is between 
~ 4.15 and 4 .87 Rq (slightly lar ger than the value of 3.65 
Rq found by IBatalha et al .1120121 ). 



9 We use the term validated to mean that the planet candidate 
is most likely a planet based on our observations but it is not a 
confirmed planet because there is not independent evidence (e.g., 
Doppler or transit timing variations) for the planet model. 



6.3 KOI 526.01 

Despite only observing a partial transit with the GTC (as il- 
lustrated in Figure[8]), we find that the transit depths clearly 
match each other, with planet-star radius ratios that are 
consistent within ~ 1.2a. While the Kepler depth is much 
shallower, there appears to be no significant contamination 
from nearby stars in the Kepler photometry (see footnote 
in H6.ip . Given that even the transit depths in the Kepler 
and the blue GTC bandpass are significantly different, and 
that dilution appears to not be the source, we believe this 
object warrants further investigation beyond the scope of 
this paper. For instance, it is possible that stellar variability 
in the target star or a blended star impacted both the GTC 
and Kepler's measured transit depths[3 Regardless, since 
the depths between the GTC bandpasses are consistent, this 
target still passes our validation test. As in the case for KOI 
420.01, we found no stars within ~ 20 arcsec of the tar- 
get that showed a transit signal at the expected time of the 
transit, so we also validate KOI 526.01 as a planet. Table [3] 
includes an updated ephemeris and impact parameter (recall 
that the transit duration was held fixed due to fitting only 
a parti al transit). Again, as for KOI 420.01, IBorucki et all 
|2011bl ) f ound an uncer t ainty of 1 on the impact param- 
eter, and iBatalha et alj l|2012T ) report an impact parame- 
ter of 0.80±0.34, so our observations and models provide a 
much stronger constraint on the impact parameter, which 
we find has an associated uncertainty of 0.0771. Based on 
our measured planet-star r adius ratios and assu ming a stel- 
lar radius of 0.92 Rq (from IBatalha et al.ll2012|. and slightly 
larger than the value of 0.80 Rq reported by Borucki et al.l 
l2011bl ). we calculate that KOI 526.01 has a radius between 
~ 6.37 and 7.16 Rq. We note that this is over tw ice the 
radius of 3.11 Rq measured by IBatalha et al.1 (|2012i ). 



6.4 KOI 1187.01 

As in the case of KOI 225.01, the light curves for KOI 
1187.01 (shown in Figure O appear to be fairly V-shaped. 
Furthermore, visual inspection of the GTC light curves 
shows that there is clearly a significant difference in the GTC 
transit depths. However, we find that the best-fit planet- 
star radius ratios only differ at a level of ~ 1.8(7, likely due 
to the large uncertainty in the radius ratio measured for 
the 858 nm light curve. We attribute this large uncertainty 



10 Kepler light curves for this target show baseline variability at 
a level of -1%. 
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to a combination of the relatively poor photometric preci- 
sions achieved for this target combined with the degener- 
acy between the impact parameter (measured to be nearly 
equal to 1) and the planet-star radius ratio. To reconcile 
these measurements with what we find visually, we com- 
pute weighted mean colors and their uncertainties for the 
in-transit and out-of-transit data. We find a mean in-transit 
color of -0.00425±0. 00050 and a mean out-of-transit color 
of -0.000160±0. 000212, which differ significantly at a level 
of ~5.8<r. Therefore, despite the consistent measured radius 
ratios, there is an obviously significant color change during 
the transit event, so we argue that KOI 1187.01 is in fact a 
false positive and not a planet. 

Contrary to KOI 225.01, the transit of KOI 1187.01 
is deeper in the redder GTC bandpass than in the bluer 
GTC bandpass. This implies that during the eclipse, more 
red light is blocked than blue light, and the smaller (eclips- 
ing) component is bluer than the primary star. In this case, 
KOI 1187 may consist of an evolved giant star that is red- 
der and severa l magn itudes brighter than the eclipsing star. 
ISlawson et al.l (|201ll ) also list this target in the eclipsing 
binary catalog as a rejected KOI (which we were also not 
aware of prior to observing t his target). Also, just as for KOI 
225.01 . ISlawson et alj (|201ll > found an orbital period for KOI 
1187.01 that is twice as long as that found bv lBorucki et al] 
(|2011bh . As discussed in Q6.1\ above and below in ^3 such 
results have ramifications on any correlation between the 
orbital period and false positive rate. 



7 DISCUSSION 

Of the four Kepler planet candidates presented in this pa- 
per, we identified two as false positives and provide further 
evidence supporting the planet hypothesis for two candi- 
dates. This suggests a false positive rate that is much higher 
than has been previously predicted for the Borucki et al. 
(2011b) KOI catalog (Morton & Johnson 2011; see TTTTIfor 
further discussion), so we consider what could cause such a 
large false positive rate for our sample. Referring back to 
Figure [T] we see that our two false positives, KOI 225.01 
and KOI 1187.01, have shorter orbital periods than the two 
KOIs we validated as planets. Therefore, our findings sug- 
gest that the false positive rate for the Kepler sample varies 
significantly with orbital period and is largest at the short- 
est periods (P < 3 days), which is what one would expect 
a priori due to the rise in the number of detached eclipsing 
binaries at these short periods (Figure [l}. 

We find no significant correlation between the false pos- 
itive rate and planet radius or other properties of the host 
star. Figure [5] illustrates that there seems to be no trend 
in the false positive rate with planet radius, as the KOIs 
we identified as false positives have the smallest and largest 
apparent radii (as measured by Kepler) in our sample. Sim- 
ilarly, there appears to be no trend with the effective tem- 
perature of the host star, although we note that both false 
positives were also the faintest targets in our sample, as il- 
lustrated in Figure [4] (also see i|7.1[) . Finally, we emphasize 
that based on the stellar parameters from the KIC, all our 
targets are likely FGK dwarfs, which are the primary targets 
of the Kepler mission. 

One factor that must be considered when making such 



conclusions is that we do not take into account the level of 
uncertainty in the stellar parameters that are drawn from 
the KIC. Any uncertainty in the stellar parameters would 
obviously affect the distribution of the planetary radii and 
stellar magnitudes and temperatures discussed here. How- 
ever, if the stellar radii in the KIC have been systematically 
over- or under-estimated, then all planetary radii would sim- 
ply scale up or down, and our conclusions about the lack of a 
correlation between the false positive rate and planet radius 
would remain the sameF^l 

Additionally, there is the issue that if a given planet 
candidate is actually an eclipsing binary with an orbital pe- 
riod that is twice as long as was initially expected, then this 
would imply that we are not probing planet candidates with 
orbital periods of less than 6 days, but instead we are ac- 
tually probing a sample of eclipsing binaries with periods 
as long as 12 days. However, since there is still a compa- 
rable population of eclipsing binaries and planet candidates 
with longer periods (e.g. out to at least 10 days, as illus- 
trated in Figure[T]), we believe this period discrepancy would 
not affect our conclusion that a population of short-period 
eclipsing binaries can significantly contaminate short-period 
planet candidates. Candidates that are initially identified as 
planets but are actually eclipsing binaries will tend to have 
shorter apparent periods (their periods will double when 
they are identified as eclipsing binaries), which strength- 
ens our argument for contamination of the planet sample 
at short periods, since this effect would further reduce the 
ratio of planets to eclipsing binaries at periods of less than 
approximately 3 days. 

Finally, we note that with the technique presented here, 
it is not possible for us to identify cases where a planet can- 
didate is actually a binary star composed of two stars that 
have the same temperature, as there would be no measur- 
able color change during the transit event. This suggests that 
some of our potentially validated planets could still be false 
positives. However, Kepler's photometry should be able to 
distinguish if the transit depths differ between every other 
transit (i.e. odd-even transit depths), which would identify 
a candidate as a false positive. In cases where a significant 
odd-even ratio is measured, this can also be used to recon- 
cile the issue with the orbital periods described above. To 
the best of our knowledge, no significant odd-even ratio was 
found for our validated planets, KOI 420.01 and KOI 526.01. 

Recently, a new list of Kepler planet candidates was an- 
nounced, bringing th e total number of candidates to 2,321 
|Batalha et alj|2012j '). As illustrated in Figures [0-131 and [1171 
the new catalog follows the same general distribution of 
orbital periods, planet radii, stellar magnitudes, t e mpera - 
tures, and Galactic latitudes as the iBorucki et al.l (|2011bT ) 
catalog, but there is notably a greater number of smaller 
and shorter-period planet candidates in the new catalog. 
This serves to emphasize the need to observationally con- 
strain the false positive rate for such small, short-period 
planets. We note that this new catalog was generated us- 
ing improved vetting metrics, so it should have a higher 



11 Some high signal-to-noise ratio eclipsing binaries are more 
likely to have correctly determined radii, so some short-period 
planet candidates that are actually eclipsing binaries are corre- 
spondingly more likely to have accurate radii. 
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fidelity than previous ones (i.e., the Borucki et a l. 2011b 
catalog). However, the KOIs from the lBorucki et all l|2011Vj ) 
catalog have not been vett e d aga inst the improved metrics 
described in lBatalha et al.l l|2012l ). so both KOI 225.01 and 
KOI 1187.01 are still included in the updated catalo g as po- 
tential planet candidates. Also. lBatalha et alj (|2012T ) remark 
that potential KOIs are not vetted against the shape of their 
light curves, so planet candidates with V-shaped transit light 
curves (such as KOI 225.01 and KOI 1187.01) are not im- 
mediately rejected. Given ongoing par allel efforts to identify 
planet candidates and KOIs and that ISlawson et al.l (|201ll ) 
previously identified both KOI 225.01 and KOI 1187.01 as 
eclipsing binaries, we recommend that the eclipsing binary 
catalog be consulted before following up any given planet 
candidate. 



7.1 Comparison to Theoretical Studies 

A recent study bv lMorton fc Johnson! ([2011) provided theo- 
retical esti mates of the false positive rate for Kepler planet 
candidates. iMorton fc Johnson! (|201 lh specifically estimated 
that nearly 90% of the 1,235 candidates presented by 
iBorucki et al.1 (|2011lJ ) had a false positi ve probability of less 
than 10%. For our targets in particular. IMorton fc Johnson! 
(|201lT ) computed the following false positive probabilities: 
0.01 (KOI 225.01), 0.04 (KOI 420.01), 0.03 (KOI 526.01) 
and 0.03 (KOI 1187.01). It is interesting to note that the 
target with the lowest false positive probability is one that 
ended up being a false positive. However, we note t h at the 
false positive probabilities from IMorton fc Johnson! (|201 lh 
are only valid for non V-shaped transit signals. Therefore, 
considering that both KOI 225.01 and KOI 1187.01 ap- 
pear to have somewhat obviously V-shaped transits, the 
lo w false positive probabiliti es computed for these targets 
bv lMorton fc Johnson (|201 lh are not appropriate. A recent 
study by IMorton! (j20ll)" improves upon the analysis from 
IMorton fc Johnson! ( 201 if ) by taking the transit shape into 
account when computing the false posit ive prob a bility for a 
given target. From their new analysis, IMorton! (|2012h cal- 
culate a false positive probability of > 0.99 for KOI 225.01 
and 0.76 for KOI 1187.01, which is consistent with our ob- 
servations. The caveat about the transit shape also affects 
the ove rall assumed false p o sitive rate for the Kepler sam- 
ple, as IMorton fc Johnson! (|201ll ) did not separate planet 
candidates with V-shaped transit signals from those with 
non- V-shaped transit signals. This means that we should 
not interpret their findings as 90-95% of Kepler planet can- 
didates are planets. However, we can use their averaged false 
positive probabilities as a starting point for the overall false 
positive rate for the Kepler sample. Despite any caveats, we 
estimate the probability of detecting 2 or 3 false positives 
out of 4 or 5 targets observed. Assuming a binomial distri- 
bution and a false positive rate of 1 0% for Kepler planet 
candidates (jMorton fc Jo hnson 201l]), we estimate that the 
probability of detecting 2 false positives from a sample of 
4 targets is less than 1% (and less than 5% for detecting 3 
false positives from 5 targets). Thus, we use this to empha- 
size that the false positive rate for current Kepler planet 
candidates with radii less than ~ 5 and orbital periods 
less than ~ 6.0 days is likely much higher than 10%. 

In regards to the full Kepler sample, 

IMorton fc Johnson! |201lf ) found that the false posi- 



tive probability varied with the depth of the transit as well 
as the magnitude and Galactic latitude of the target star, 
but they did not investigate how the probability might 
vary with orbital period or transit duration (though they 
claim that such properties would serve to decrease their 
probability estimates and thus their estimates are upper 
limits). Specifically, our results support their finding that 
the false positive rate increases slightly with magnitude, as 
our two faintest targets ended up being false positives. They 
also find a general increase in the false positive probability 
with increasing transit depth (though there are some 
local minima and maxima; e.g., see Figure 7 in Morton & 
Johnson 2011) and with decreasing Galactic latitude. In 
our case, the target with the largest transit depth (as mea- 
sured in the Kepler bandpass) was KOI 420.01, which we 
validated as a planet. While we find no obvious correlation 
in our sample between Kepler's measured transit depths 
and the false positive rate, dilution is an important factor 
that has to be considered (e.g., as illustrated by our light 
curves for KOI 225.01 in Figure |5J, as some of the depths 
measured by Kepler may be underestimated. In regards 
to Galactic latitude, in Figure \W\ we show t he cumulative 
distribu tion of the Gal actic latitudes f or th e IBorucki et al.l 
ll2011bl) KOI list t he iBatalha et all (|2012l ) KOI list and 
the Slawson et al.l (|2011 ) catalog of eclipsing binaries, 
which illustrates that the samples have the same general 
distribution. Furthermore, our false positives are indeed 
located at lower lat i tudes , which supports the argument by 
IMorton fc Johnson! l|201ll) that there is a slight increase in 
the false positive rate with decreasing Galactic latitude. 
We present a similar figure in Figure 111! where we show 
the cumulative distribution of the Galactic latitudes for the 
Kepler eclipsing binaries along with separate distributions 
for stars that host short-period (P ^ 2 days) and long- 
period (P > 2 days) planet candidates. We see that the 
distribution for the stars that host short-period candidates 
matches the eclipsing binary distributions more so than 
the sample of long-period candidates. Therefore, Figure [TT] 
further illustrates the likely contamination of short-period 
planet candidates by eclipsing binaries. In Figure 1121 we 
illustrate how the number of eclipsing binaries varies with 
Galactic latitude, and we find that overall the fraction of 
Kepler targets that are eclipsing binaries is about the same 
as a function of Galactic latitude. This further implies a 
nearly consistent presence of eclipsing binaries that could 
plague the Kepler planet candidate sample. 

7.2 Comparison to Observational Studies 

IDesert et all (|2012l ) have used waim-Spitzer follow-up of 
KOIs to find a low false positive rate c onsistent with the es- 
timates from IMorton fc Johnson! (|201lf ). Overall, our sample 
includes shorter period planets than their sample, as well as 
fainter targets (excluding their sample of M dwarf stars). 
Given that only planet candidates with the shortest orbital 
periods (<< 2 days) are likely to be significantly contami- 
nated by non-d etached eclip s ing b inaries, it is not entirely 
surprising that IDesert et alj (|2012h found a low false posi- 
tive rate for their sample of planet candidates with periods 

of ~1. 8-100 days. 

Recently, ISanterne et al.l |2012i ') presented a Doppler 
study of short-period giant Kepler planet candidates. They 
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Figure 10. Cumula tive fractio n of planet-hosting stars from 
iBorucki etldl ll2011bl , blue) and IBatalha et~all ||2012| , gray) and 
detached (red) and all "other" (primari ly semi-detached and over- 
contact; black) eclipsing binaries from lSlawson et alj d201lh as a 
function of Galactic latitude. The cumulative functions w ere nor- 
malized against 99 7 stars for the IBoru cki et al. | ll2011bh sample, 
1,790 stars for the IBatalha et al sample (which includes 

the Borucki et al. 2011b sample), 1,266 detached eclipsing bi- 
naries, and 901 "other" eclipsing binaries (from Slawson et al. 
2011). The gray vertical dashed lines indicate the Galactic lati- 
tude of the targets we observed with the GTC, with each target 
marked accordingly. The rate of planet-hosting stars is slightly 
larger towards higher latitudes, while the rate of eclipsing bina- 
ries is slightly larger at lower latitudes, but in general the two 
samples are not significantly different. Both the KOIs we found 
to be false positives are located at lower latitudes. 

targeted a sample of 46 Kepler planet candidates with 
transit depths greater than 0.4%, orbital periods less than 
25 days, and host stars brighter than 14.7 magnitude 
(Kp). Based on the i r radi al velocity follow-up observa- 
tions, ISanterne et alj (|2012l ) estimate a false positive rate 
of 34.8±6.5% for this subset of Kepler planet candidates. 
While we probed a different (and smaller) population of 
Kepler planet candidates, our observed false positive rate 
of 50% (with large uncertainty due to our small sample 
size) for small, short-period Kepler planet ca ndidates is sup- 
ported by their observations. The study by ISanterne et al.1 
|2012r ) also supports our argument that different popula- 
tions of Kepler targets likely have different false positive 
rates associated with them. 



8 CONCLUSION 

We acquired multi-color transit photometry of four small 
(~2.5-4.9 R ffi ), short-period (~0.37-6.0 days) Kepler 
planet candidates with the GTC. Based on the transit color, 
we identified two candidates as false positives. For two, we 
find further evidence supporting the planet hypothesis, con- 




Galactic Latitude [degrees] 

Figure 11. Similar to Figure [lOl Here, we show the cumulative 
fraction of stars hosting short-period (P ^ 2 days; long-dashed 
blue curve) and long-period (P > 2 days; dotted blue curve) KOIs 
(based on the full list of 2,321 KOIs from Batalha et. al 2012) 
along with the cumulative fract ion of two eclipsing binary popula- 
tions from lSlawson et al. shown in Figure IT7J1 as a function 
of Galactic latitude. The curves for the planet-hosting stars were 
normalized against stars hosting 158 short-period and 2,163 long- 
period planet candidates. We find that the short-period planet 
population is more consistent with the eclipsing binary popula- 
tions than the long-period planet population, which further sup- 
ports our argument for contamination of the planet population 
by eclipsing binaries at short periods. Conversely, the long-period 
planet population is less likely to be infiltrated by eclipsing bina- 
ries. See the caption to Figure [Tol for additional details. 



sistent with validated planets. We also remind the reader 
of KOI 565.01, which was a planet candidate observed in a 
similar fashion and that was also found to be a false positive 
(albeit it was selected from the first KOI catalog published 
by Borucki et al. 2011a; Colon & Ford 2011). While we find 
a high false positive rate (2/4 or 3/5, if we include KOI 
565) in our small sample, we caution that this is likely not 
representative of the entire sample of Kepler planet candi- 
dates, due to the small number of targets we observed and 
the specific properties of these candidates (e.g. the orbital 
period and size). Nevertheless, our results demonstrate the 
importance of considering these properties when evaluating 
the false positive probability of specific systems. While our 
findings seem to contradi ct the theoretical estimates from 
iMorton Sz Johnson! (|201ll ). the low false positive rate that 
they estimate is based on the assumption that all candi- 
dates had passed preliminary false-positive vetting metrics 
based on Kepler photometry and astrometry. Thus, if we 
consider the obviously V-shaped transits for KOI 225.01 and 
KOI 1187.01 to imp ly a false positive nature f or these KOIs, 
then according to IMorton &; Johnson! (|201ll ) the low false 
positive probabilities for these targets are not accurate. The 
false positive rate for our sample is also much larger than 
the observational constraints from iDesert et all |2012f ) that 
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Figure 12. Histograms of the numbe r of detached ( r ed) an d all 
"other" (gray) eclipsing binaries from IsTawson et al.l (|201lh as a 
function of orbital period. Each panel shows a different bin of 
Galactic latitude. As in Figure [l] the righthand y-axis shows the 
percentage of systems relative to the total number of target stars 
observed by Kepler in Ql. There is a roughly consistent number 
of eclipsing binaries in each latitude bin (except for the highest 
latitudes), which suggests a consistent presence of eclipsing bi- 
naries that could contaminate potential planet candidates. While 
there are fewer eclipsing binaries at the highest latitudes, there 
are also correspondingly fewer planet candidates at those lati- 
tudes (see Figure [lO}. 



predict that the false positive rate is much less than 10%. 
This is likely partly a result of different targe ts being probed 
by the different studies. The recent study bv lSanterne et all 
l|2012l ) further supports this idea, as they found a ~35% 
false positive rate for short-period giant Kepler planet can- 
didates. We plan to continue observing small, short-period 
KOIs with the GTC in order to improve the sample size of 
our study. The observations presented here, as well as future 
observations with the GTC, greatly complement follow-up of 
KOIs done with warm- Spitzer as well as other observatories. 
As we can expect the number of Kepler planet candidates to 
continue to increase, we can use results from all such studies 
to pinpoint which targets are the best to follow-up in order 
to maximize the science output from the Kepler mission. 



los Muchachos of the Instituto de Astroffsica de Canarias, 
on the island of La Palma. The GTC is a joint initiative of 
Spain (led by the Instituto de Astroflsica de Canarias), the 
University of Florida and Mexico, including the Instituto 
de Astronomia de la Universidad Nacional Autonoma de 
Mexico (IA-UNAM) and Instituto Nacional de Astroflsica, 
Optica y Electronica (INAOE) . 
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